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It is shown that bifurcations of the mean-field dynamics of a Bose-Einstein condensate can be 
related with the quantum phase transitions of the original many-body system. As an example we 
explore the intra-band tunneling in the two-dimensional optical lattice. Such a system allows for 
easy control by the lattice depth as well as for macroscopic visualization of the phase transition. 
The system manifests switching between two selftrapping states or from a selftrapping state to 
a superposition of the macroscopically populated selftrapping states with the step-like variation 
of the control parameter about the bifurcation point. We have also observed the magnification 
of the microscopic difference between the even and odd number of atoms to a macroscopically 
distinguishable dynamics of the system. 
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Introduction.- Since the very beginning of the quan- 
tum mechanics its relation to the classical dynamics con- 
stitutes one of the central questions of the theory. De- 
pendence of the energy levels distribution on the type 
of dynamics of the corresponding classical system [l[ , in 
general, and the quantum system response to variation 
of the bifurcation parameters controlling the qualitative 
changes of the classical behavior 0] are among the ma- 
jor issues 0. One of the main tools in studies of the 
quantum-classical correspondence is the WKB approxi- 
mation, where, loosely speaking, the Planck constant h 
is regarded as a small parameter. 

On the other hand, for a iV-boson system the limit 
N — > oo at a constant density, leading to the mean-field 
approximation, can also be understood as a semiclassi- 
cal limit. This latter approach has received a great deal 
of attention during the last decade Q , due its high rele- 
vance to the theory of Bose-Einstein condensates (BECs), 
many properties of whose dynamics are remarkably well 
described within the framework of the mean-field mod- 
els [B[ . More recently, it was shown [f| 0] that the mean- 
field description of a few-mode iV-boson system can be 
recast in a form similar to the WKB approximation for a 
discrete Schrodinger equation Q, emergent for the coef- 
ficients of the wavefunction expansion in the associated 
Fock space, where 1/N plays the role similar to that of 
the Planck constant in the conventional WKB approxi- 
mation. 

The mean-field equations of a system of interacting 
bosons are nonlinear, hence, they naturally manifest 
many common features of the nonlinear dynamics, in- 
cluding bifurcations of the stationary solutions caused by 
variation of the system parameters. One of the well stud- 
ied examples is a boson- Josephson junction [§], which can 
show either equally populated (symmetric) or strongly 



asymmetric states, characterized by population of only 
one of the sites (the well known phenomenon of selftrap- 
ping Now, exploring parallels between the semi- 
classical approach and the mean-field approximation one 
can pose the natural question: what changes occur in a 
manybody system when a control parameter crosses an 
instability (e.g. bifurcation) point of the limiting mean- 
field system? 

In the present Letter we give a partial answer showing 
that one of the possible scenarios is the quantum phase 
transition of the second type, associated with the switch- 
ing of the wave-function in the Fock space between the 
"coherent" and "Bogoliubov" states possessing distinct 
features. Considering a flexible (time-dependent) control 
parameter, we have also found a strong sensitivity of the 
system to the parity of the total number of atoms N, 
showing parity-dependent structure of the energy levels 
and the macroscopically different dynamics for different 
parity of N. Observation of the discussed phenomena is 
feasible in the experimental setting available nowadays. 

Quantum and mean-field models.- We consider the 
nonlinearity-induced intra-band tunneling of BEC be- 
tween the two high-symmetry X-points of a two- 
dimensional square optical lattice (OL). The process is 
described by the two mode boson Hamiltonian (see 
for the details) 
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where 6 7 - and b\ are the annihilation and creation oper- 
ators of the two X-states, A (0 < A < I) is the lat- 
tice parameter easily controllable by variation of the lat- 
tice depth (or period). The Schrodinger equation for 
the BEC in a state |*) reads ihd T \%) = H\%), where 
h = 2/N and t = (2gp/h)t, with g = 4-Kh 2 a s /m 
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and the atomic density p. The link with the semiclas- 
sical limit is evident for the Hamiltonian in the form 
(fTJ): the Schrodingcr equation written in the Fock ba- 

(b\) k (bt) N ~ k 
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:|0), depends only on the 
- k)/N, while h serves 
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relative populations k /N and (N 
as an effective "Planck constant" . 

Hamiltonian |T]) represents a nonlinear version of the 
well-known boson- Josephson model (see, e.g. 0, [111]), 
where unlike in the previously studied models the states 
are coupled by the exchange of pairs of atoms. This is 
a fairly common situation for systems with four-wave- 
mixing, provided by the two-body interactions involving 
four bosons. The exchange of the bosons by pairs results 
in the coupling of the states with the same parity of the 
population and is reflected in the double degeneracy of all 
(N + l)/2 energy levels for odd N, due to the symmetry 
relation (2k, N - 2k\^i) = (N - 2k, 2fc|* 2 )< For even N 
the energy levels show quasi degeneracy (see below). 

The mean-field limit of the system ([1]) can be formally 
obtained by replacing the boson operators bj in ([1]) by the 
c-numbers bi — > \J Nxe l ^/ A and b 2 
what gives the classical Hamiltonian 
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H = x(l- x) [2A- 1 + Acos( 
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where x = (n-y)/N is the population density and (j> = 
arg((&2) 2 &i) is the relative phase. H possesses two sta- 
tionary points describing equally populated AT-statcs: 
the classical energy maximum Pi = (x = \ , 4> — Q) and 
minimum Pi = (x = h , <p = tt). Pi is dynamically stable 
in the domain A > A c = |. For A < A c it looses its 
stability, and another set of stationary points x = 1 (Si) 
and x = (S2) appears, which is a fairly general situa- 
tion in nonlinear boson models. The appearing solutions 
describe the symmetry breaking leading to selftrapping. 

Energy levels near the critical point.- To describe the 
spectrum of the Hamiltonian JTJ) in the vicinity of the 
critical value A c we rewrite H in terms of the operators 
01,2 = {h Tih)/V2 

H = H +(l-^)v+£(A), £(A) = ±±±+± (3) 



where Hq = j^a\aia\a2 and V — (a\ci2 + a 2 a\ 
At the critical point the energy spectrum is determined 
by H : E m = jj£m(N - m) + £ (A c ), where m is the 
occupation number corresponding to the operator a\a\. 
The spectrum of Hq is doubly degenerate (except for the 
top level for even N) due to the symmetry m — > N — m. 
The ground state energy is £ m i n (A c ) = Eq = En, while 
the top energy level has m — N/2 for even N and 
to = (N ± l)/2 for odd N. Restricting ourselves to even 
number of bosons we get £ m ax{A c ) = \ + 4w- 



Now consider small deviations of A from the bifurca- 
tion point A c . To this end, for a fixed N, one can use 
the basis consisting of the degenerate eigenstates of H : 

(<4) m (4 i) N ~ m 

■ 0), j = 1,2, m = 0,..,f. 
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The conditions for V to be treated as a perturbation de- 
pend on to as is seen from the diagonal matrix elements: 



(E m ,j\V\E m ,j) 
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At the lower levels (m <C N/2) the energy gaps between 
the degenerate subspaces and the perturbation both scale 
as AP ~ ^V _1 , hence the condition of applicability is 
I A — A c I -C 1 and the lower energy subspaces acquire 
simple shifts. At the upper energy levels (m ~ N/2) the 
above energy gaps behave as AP ~ Af -2 . Since (V) ~ 1 
in this case, the perturbation theory is applicable only in 
an interval of A of the size on the order of A^~ 2 . There is 
a dramatic transition in the energy levels, e.g. Fig. [TJb) 
shows the exchange of the double degeneracy of the top 
levels for even N in this Af~ 2 -small interval of A. By 
considering the phase of 



(E m ,j\(blbif\E m ,j)=- 
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-to (N — to) , (5) 



it is easy to verify that the upper and lower eigenstates 
correspond, respectively, to the mean-field stationary 
points Pi (<f) = 0) and P 2 (</> = 7r). 
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FIG. 1: (a) The energy levels of H for N = 200 and (b) a de- 
tailed picture in the vicinity of Ac. The classical energy lines 
of the mean-field fixed points Pi and & are visibly formed. 
The top energy levels for sufficiently large |A — A c | are quasi- 
degenerate with the inter-level distances indistinguishable on 
the scale of the figure (see the discussion in the text below). 



Spectrum in the limit N — > oo. Coherent states and 
selftrapping states. - For A^ 1 — A -1 3> A^~ 2 the quantum 
states corresponding to Pi can be obtained by quantizing 
the local classical Hamiltonian ([2|), i.e. by expanding it 
with respect to x— 1/2 and <fr and setting <j> = —ih-J^ (see 
also Ref. [l2| ; on this way one looses the term of order 
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1/N in £ m ax)- The "wave function" tp(x) 
VN (k, N - k\ip) satisfies 
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ip = Eip. (6) 



Eq. © is the negative mass quantum oscillator prob- 
lem with the frequency to 2 = 8(3— -^). The respec- 

p(top) 



tive descending energy levels read E\ ' = £ max + 
i (4- - - ^ (n + |) . The eigenfunctions are local- 

eigenfunction is 
In the original discrete 
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ized in the Fock space, e.g. the n 
tp (x) = Cexp (x - \) 2 
variable x = k/N, there are even and odd eigenstates 
C2k and Cik-\ related by the approximate symmetry 
C n ~ Cn+i, hence the energy levels are quasi doubly 
degenerate [c.f. Fig. 1(b)]. 

The local approximation becomes invalid as 
A^ 1 — A -1 ~ A^~ 2 (the wave-function delocalizes). 
The other set of the stationary points, Si, 2, becomes 
stable for A < A c in the mean-field limit. In this case, 
however, the phase <f> is undefined. Let us first consider 
the full quantum case, for example, the limit (ni) <C N 
(i.e. the point S2). The resulting reduced Hamiltonian 
can be either easily derived in the Fock basis or ob- 
tained by formally setting 62 = N and retaining the 
lowest-order terms in bi and b\: 



H ~ H S2 — - 
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Hamiltonian {7} can be diagonalized by the Bogoliubov 
transformation c = cosh(6>)&i — smh(6)b^ , where 9 = 
9(A) > is determined from tanh(20) = A/(l - 2A). We 
get 



Hs 2 — - 
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Thus c^c gives the number of negative-energy quasi- 
particles over the Bogoliubov (squeezed) vacuum solv- 
ing c|vac) = 0. In the atom-number basis |vac) is 
a superposition of the Fock states with = 

0, (Cq is a nor- 



tanh fc (0)7(2l)!/(2 fc fc!)Co and 



malization constant). 

The validity condition of the approximation [[7]). 
given by (n),An <C N, can be rewritten in the form 
tanh _2 (26*) > 1 + A^~ 2 , what is the same as A^-A^ 1 > 
A^ -2 . In this case, the eigenstates of (jTJ) are well- localized 
in the atom-number Fock space, i.e. the coefficients Cik 
decay fast enough. The condition for this excludes the 
same small interval as in the perturbation theory, hence 
the transition between the coherent states and the self- 
trapping (Bogoliubov) states occurs on the interval of 
A of order of A^~ 2 . The convergence of the eigenstates 
of Hs 2 to that of the full Hamiltonian JT]) turns out 
to be remarkably fast as it is shown in Fig. [DJa). In 



Fig. [2}d the dramatic deformation of the top energy eigen- 
state of H (corresponding to the S2-P1 transition) about 
the critical A c is shown. Finally, we note that for even 
N the quasi double degeneracy of the energy levels for 
A^-A^ 1 > N~ 2 (c.f. Fig. 1(b)) is due to the exchange 
symmetry between Si and S2 resulting in equal energy 
levels of the Hamiltonians Hs 1 and Hs 2 ■ 



0.338 



0.336 



0.332 





10 15 20 25 30 35 40 
k 



FIG. 2: (Color online) (a) Convergence of the four upper 
eigenstates of the Hamiltonian |T} to the eigenstates of Hs 2 
(shown by dots) for for N = 100 (solid lines) and N = 1000 
(dashed lines), for A = A c — 0.1. (b) The contour plot of the 
state corresponding to the top energy level in the vicinity of 
A c for N — 200. 



In the mean-field description of the stationary point 
S2 the associated Hamiltonian is defined by replacing the 
boson operators in Eq. (JTJ by the c- numbers b\ = \/~N a 
and 62 = y/N/3. Using \a\ 2 + \(3\ 2 = 1 and fixing the 
irrelevant common phase by setting f3 real we get the dy- 
namical variables a and a* and the classical Hamiltonian 
in the form Hs 2 = \ + \{l ~ |a| 2 ){2(2A - l)|a| 2 + A[a 2 + 
(a*) 2 ]}, from which the stability of the point S2 (a = 0) 
for A < A c follows. 

Thus, the passage through the bifurcation point A c of 
the mean-field model, corresponds to the phase transition 
in the quantum many-body system on an interval of the 
control parameter scaling as N~ 2 and reflected in the de- 
formation of the spectrum and dramatic change of the 
system wave-function in the Fock space. The described 
change of the system is related to the change of the sym- 
metry of the atomic distribution, and thus it is the second 
order phase transition. 

In our case this scenario corresponds to loss of stability 
of the selftrapping solutions Si and S2 and appearance of 
the stable stationary point Pi. In the quantum descrip- 
tion this happens by a set of avoided crossings of the top 
energy levels (and splitting of the quasi-degenerate en- 
ergy levels for even N) as the parameter A sweeps the 
small interval on the order of A^~ 2 about the critical 
value A c (see Fig.[T]). For lower energy levels the avoided 
crossings appear along the two straight lines approximat- 
ing the classical energies of the two involved stationary 



points: H{Pi) = \ + 
(A > Ac), see Fig. [J 
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Dynamics of the phase transition.- Let us see how the 
quantum phase transition shows up in the system dynam- 
ics when A is time-dependent. The selftrapping states Si 
and 52, eigenstates of the Hamiltonian fl}, correspond 
to occupation of just one of the A-points. Such an ini- 
tial condition can be experimentally created by switching 
on a moving lattice with A < A c (see e.g. [7|). As the 
lattice parameter A(r) passes the critical value from be- 
low, the selftrapping states are replaced by the coherent 
states with comparable average occupations of the two 
A-points. 

A more intriguing dynamics is observed when A(r) is 
a smooth step-like function between Ai and A2 such that 
Ai < A c < A2. In this case, the system dynamics and 
the emerging states dramatically depend also on parity 
of the number of atoms. For fixed Ai^ the system be- 
havior crucially depends on the time that A(r) spends 
above A c . More specifically, one can identify two dis- 
tinct scenarios, which can be described as a switching 
dynamics between the selftrapping states at the two Ap- 
points, Fig. [U(a),(b) or dynamic creation of the super- 
position of macroscopically distinct states, well approx- 
imated by {Ck\k, N- k) + C N - k \N - k,k)) with a 
k<k m 

small k m /N, Fig. 0)[c), (d) (where k m /N m 0.2). In the 
case of macroscopic superposition the dynamics shows 
anomalous dependence on parity of N , i.e. showing the 
same behavior for large N of the same parity but macro- 
scopically distinct behavior for N and N + 1, Fig. [3Jc). 
Note that the mean-field dynamics is close to the quan- 
tum one in the switching case, Fig.[3la), while it is dra- 
matically different in the superposition case, Fig. [3jc) . 




50 100 150 200 50 100 150 200 




50 100 150 200 50 100 150 200 

i x 



FIG. 3: (Color online) The average population densi- 
ties {m)/N, (a) and (c), and the atom- number prob- 
abilities \C k \ 2 , (b) and (d), for A(r) = Ai + (Aa - 
Ai) [tanh (r — n) — tanh (r — m)] /2. The corresponding 
classical dynamics is shown by the dash-dot lines in (a) and 
(c). Here Ai = 0.25 < A cr , A 2 = 0.5 > A cr , t\ = 50 and 
T2 = 85 (a) with N = 500 and 501 (indistinguishable), while 
in (c) T2 = 135 with TV = 500 and 400 (the upper solid and 
dashed lines) and N = 501 and 401 (the lower lines). The 
initial state is |vac) of Hs ± , but using |AT, 0) gives a similar 
picture. 



To estimate the physical time scale, t = i p hT = 
5 "l d l tf t, we assume that a condensate of 87 Rb atoms 
is loaded in a square lattice with the mean density of 
7Vp C = 20 atoms per cite. If the lattice constant d — 2 //m 
and the oscillator length of the tight transverse trap (to 
assure the two-dimensional approximation) £j_ — 0.1 /xm, 
then t p h ~ 0.2 ms and the time necessary for the creation 
of the macroscopic superposition of Fig. Etc),(d) is about 
20 ms. 

Conclusion.- We have shown that behind the mean- 
field instability in the intra-band tunneling of BEC in 
an optical lattice is a quantum phase transition be- 
tween macroscopically distinct states, giving a macro- 
scopic magnification of the microscopic quantum features 
of the system. A spectacular demonstration of this is the 
dynamic formation of the superposition of macroscopi- 
cally distinct states, which, besides being responsible for 
the difference between the mean-field and quantum dy- 
namics (see also recent Ref. [13|]), shows also an anoma- 
lous dependence on parity of BEC atoms reflecting dis- 
tinct energy level structure for even and odd number of 
atoms. 
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